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jj 要 : 差分 进化 是 一 种 有 效 的 优化 技术 ， 已 成 功用 于 多 目标 优化 问题 。 但 也 存在 Pareto 最 优 集合 的 收敛 慢 和 多 样 性 


差 等 问题 。 针 对 上 述 不 足 ， 提 出 了 一 种 基于 分 解 和 多 策略 变异 的 多 目标 差分 进化 算法 (MODE/DMSM)。 该 算法 利用 基 
于 分 解 的 方法 将 多 目标 优化 问题 分 解 为 多 个 单 目 标 优化 问题 通过 高 效 的 非 支 配 排 序 方法 选择 具有 良好 收 伊 性 和 多 样 
性 的 解 来 指导 差分 进化 过 程 ; 采用 了 多 策略 变异 方法 来 平衡 进化 过 程 中 收敛 性 和 多 样 性 。 在 ZDT 和 DTLZ 的 10 个 测 


试 函 数 上 的 仿真 结果 表明 ， 本 文 算 法 在 Parato 最 优 集合 的 收敛 性 和 多 样 性 优 于 其 他 六 种 代表 性 多 目标 优化 算法 。 
关键 词 : 多 目标 优化 ; 差分 进化 ; 分 解 ; 多 策略 变异 
中 图 分 类 号 : TP301.6 doi: 10.3969/j.issn.1001-3695.2018.01.0028 


Multiobjective differential evolution algorithm based on 
decomposition and multi-strategy mutation 
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Abstract: Differential evolution algorithm is an efficient optimization technique that has been successfully applied to 
multiobjective optimization problems. However, there are also some defects, i. e. the slow convergence and poor diversity of 
the Pareto optimal set. Addressing these issues, this paper presented a multiobjective differential evolution algorithm based on 
decomposition and multi-strategy mutation (MODE/DMSM) . MODE/DMSM utilized the decomposition-based approach to 
decompose a multiobjective optimization problem into multiple single-objective optimization problems. Moreover, 
MODE/DMSM adopted the efficient non-dominated sorting approach to select solutions which have both good convergence and 
diversity to guide the differential evolutionary process. Eventually, MODE/DMSM employed the multi-strategy mutation 
approach to balance the convergence and diversity in the evolutionary process. The results of simulations on 10 test functions 
of ZDT and DTLZ show that MODE/DMSM outperforms than the other six representative multiobjective optimization 
algorithms in terms of the good convergence and diversity of the Pareto optimal set. 
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m 解决 多 目标 问题 。 当 前 多 目标 进化 算法 主要 分 为 三 类 : 
a) 基于 Pareto 占 优 的 多 目标 进化 算法 采用 Pareto 支配 来 
在 优化 问题 中 存在 同时 需要 优化 两 个 或 者 多 个 互相 冲突 的 ” ”选择 非 支 配 解 ， 如 NSGAII 外 和 SPEA2P!, 
标 问 题 ， 这 类 问题 称 为 多 目标 优化 问题 [Hmultiobjective b) 基 于 评估 指标 的 多 目标 进化 算法 采用 某 种 评估 指标 来 选 
optimization problem, MOP)。 多 目标 问题 没有 单个 最 优 解 , 其 最 PEdESCHOAT, Ju IBEAUTRI HypED。 
优 解 是 一 个 Pareto 解 集 。Pareto 解 集 的 收敛 性 和 多 样 性 是 衡量 c) 基 于 分 解 技术 的 多 目标 进化 算法 将 多 目标 问题 分 解 成 一 
多 目标 算法 的 重要 指标 。 组 单 目 标 优化 问题 。 如 Zhang 和 Li 提出 了 一 种 基于 分 解 的 多 
多 目标 进化 算法 (multiobjective evolutionary algorithms， 标 进 化 (MOEA/D) 算 法 器, 该 算法 将 多 目标 问题 分 解 成 多 个 单 
MOEA) 是 基于 种 群 模拟 生物 进化 全 局 优化 算法 , 已 经 广泛 应 用 目标 问题 来 求解 ， 然 后 用 基因 遗传 算法 (genetic algorithm, GA) 
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计算 


求解 这 些 单 目标 问题 。 

差分 进化 算法 中 (differential evolution; DE) 是 一 种 简单 高 
效 的 全 局 优化 进化 算法 Cevolutionary algorithm, EA), 已 经 应 
用 于 多 目标 问题 求解 。Zhong 和 Zhang 提出 了 一 种 带 有 随机 编 
码 策略 的 自 适 应 多 目标 差分 进化 算法 (AS-MODE)II， 使 用 DE 
中 一 种 简单 的 变异 策略 来 生成 新 种 群 。Venske 提出 一 种 自 适应 
的 多 目标 差分 进化 算法 (ADEMO/D)91。Jiang 提出 了 一 种 基于 两 
阶段 策略 和 小 生境 指导 方案 的 改进 MOEA/D(MOEA/D-TPN)W 9 
该 算法 中 DE 作为 基础 优化 算法 ,文献 [11] 提 出 了 基于 分 解 的 差 
分 进化 算法 (MOEADDE)W 1， 使 用 DE 代替 MOEA/D 中 GA 算 
法 ， 通 过 分 解 的 方法 更 新 种 群 。 总 体 来 看 ， 多 目标 进化 算法 主 
要 采用 单一 的 变异 策略 ， 导 致 Pareto 解 集 的 收敛 性 和 多 样 性 不 
够 理想 。 

为 了 提高 Pareto 解 集 的 收敛 性 和 多 样 性 ， 本 文 在 
MOEADDE 的 基础 上 提出 了 一 种 基于 分 解 和 多 策略 变异 的 多 


I 


目标 差分 进化 算法 (MODE/DMSM)。 算 法 采用 改进 Tchebycheff 


Z 分解 方法 将 一 个 多 目标 问题 分 解 为 一 组 单 目 标 优化 子 问题 ， 同 
O 时 生成 一 系列 均匀 分 布 的 权重 向 量 。 为 了 兼顾 收敛 性 和 多 样 性 ， 
N 算法 采用 高 效 非 支配 排序 对 种 群 进行 排序 并 保存 了 非 支 配 解 和 
CÓ 一 部 分 较 好 的 解 。 算 法 利用 多 策略 变异 差分 进化 算法 来 引导 种 
er 群 向 真实 的 Pareto 前 沿 收敛 同时 均匀 分 布 ， 实 现 收 敛 性 和 多 样 
CO 性 的 平衡 。 与 其 他 典型 的 多 目标 优化 算法 进行 仿真 比较 ， 结 果 
= 表明 本 文 提出 的 算法 在 解决 多 目标 问题 时 在 收敛 性 和 多 样 性 具 
C 有 一 定 的 优势 。 

N 

S 1 ， 相关 背景 

>< 11 多 目标 优化 问题 

E 一 般 的 多 目标 问题 , 以 最 小 化 形式 为 例 , 具有 n 个 决策 变 
-三 。 量 和 个 优化 目标 ， 可 表示 为 如 下 形式 ， 

© minimize F(x) = Cfi (x)..... fa C) 


a) 


subject to: xeQ 


1.2 基于 Tchebycheff 多 目标 分 解 方法 

分 解 方法 将 多 目标 优化 问题 转换 成 单 目标 优化 问题 , 给 每 
个 个 体 分 配 一 个 不 同 的 权重 向 量 来 缓解 目标 之 间 的 冲突 。 常 见 
的 分 解 方法 有 切 比 雪夫 (Tcehebyche 提 方法 四。 传统 的 
Tchebycheff 分 解 方法 公式 如 下 : 


T 


minimize : g(x | 4, z") = max(4, fi -=z 


! (5) 


subject to : x eQ 
其 中 : 4 为 第 i 个 权重 向 量 ; 
点 ; z-min(fo)|xeQO), 
139 差分 进化 算法 
差分 进化 算法 是 Storn 和 Price 提出 的 一 种 基于 种 群 的 随 


z = (zi,…z) 为 更 新 过 程 中 的 理想 


机 搜索 算法 ， 主 要 通过 变异 、 交 叉 和 选择 等 操作 来 进行 智能 搜 
Ro DE 与 其 他 进化 算法 的 主要 区 别 在 于 变异 操作 , 目前 被 广泛 
采用 的 变异 操作 如 下 : 
DE/rand/103]; 

vo =i, o +F: (ipe Ao) (6) 
DE/rand-to-best/103): 

Vo m3 th Ce o) FG s 7) (7) 


DE/best/103): 

"EPIIT HERO (8) 
其 中 : n, nAn EARR NP 中 随机 产生 的 3 个 互 不 相 
等 的 整数 ，G 为 当前 种 群 的 进化 代数 ，F 为 缩放 因子 ，best K 
示 当 前 最 好 个 体 的 下 标 。v 和 + 分别 表示 变异 向 量 和 种 群 的 个 
体 。 


2 ”基于 分 解 和 多 策略 变异 的 多 目标 差分 进化 算法 


2.1 多 目标 分 解 方法 
Tchebycheff 函数 是 在 基于 分 解 的 MOEAs 中 使 用 的 最 常 
见 类 型 的 聚合 函数 之 一 。 在 本 文中 采用 了 文献 [14] 中 改进 的 


Tchebycheff 函数 04， 


这 种 改进 的 Tchebycheff 函数 在 原始 函数 


上 有 两 个 优点 。 首先 


,均匀 分 布 的 权 向 量 导致 目标 空间 中 搜索 


方向 的 均匀 分 布 ; 


其 次 ， 每 个 权重 向 量 对 应 了 


F 位 于 PF 上 唯 


Ropa x-e(x.x)jeQcR" 称 为 决策 向 量 ; 
FG) - G0... CD) eR 称 为 目标 向 量 ， 为 决策 空间 。 以下。 的 解 "9。 这 两 个 优点 在 一 定 程 度 上 缓解 了 多 样 性 保存 的 困难 。 
是 多 目标 问题 引出 的 几 个 相关 定义 和 术语 。 EE Med UE 

定义 1 Pareto ZEN, 假设 ,x eog amag Minimizes D = mes en [A1 (9) 


subject to : xeQ 
其 中 : z 为 理想 点 ， 妈 为 权重 向 量 。 勾 >0， 如 果 和 =0， 将 设 
| 4 - 10^ , Das 和 Dennis0g 提 出 标准 边界 交叉 (NBI) 技术 被 
泛 的 应 用 在 目标 空间 生成 权 习 


Hx * bx Æ Pareto 支配 的 ， 记 为 过 <x， 当 且 仅 当 
Vi=1,.. m, f(x) € fŒ) ^ F) + F(x) 2) 


定义 2 Pareto ERN, Pareto 最 优 解 集 (Pareto Set, 
简称 PS) 是 所 有 Pareto 最 优 解 的 集合 。 


行 解 , 


pe rr 
可 
ji 


depo 


; ,0 H M Qj (10) 
A; dir =1 


PS ={x |23xeQ,x 4x) (3) 


3EX.3 Pareto ERTA. Pareto 最 优 前 沿 (Pareto front, 


其 中 : 将 每 维 目 标 分 割 成 五 个 相等 的 区 间 ，M 为 目标 个 数 ， 则 
PF) 是 Pareto 最 优 解 集 PS 中 的 解 对 应 目标 函数 值 组 成 的 集合 。 A I 
生成 的 权重 向 量 的 数量 N = Chima 。 假 设 将 三 维 空间 每 维 目标 

PF -(F(x )|x € PS} (4) " no utere T " 
分 割 成 2 个 相等 的 区 间 ， 生 成 权重 向 量 的 个 数 N = C6, B 


1 是 在 三 维 空间 上 生成 权重 向 量 的 示意 图 。 
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图 1 权重 向 量 生成 示意 图 
本 文采 用 了 上 述 的 方法 4 
权重 向 量 使 得 在 目 


增加 了 种 群 的 多 样 性 和 收敛 性 
2.2 ”高效 非 支配 排序 策略 


本 文 运用 了 高 效 非 支配 排序 3 来 保存 比较 优 的 解 。 这 种 方 
法 首先 根据 第 一 个 目标 值 对 种 群 己 中 的 N 个 解 进行 升序 排序 ， 


其 中 六 是 种 群 大 小 。 如 果 两 个 解 的 第 一 个 


按 第 二 个 


排序 。 如 果 解 在 所 有 目标 函数 值 中 都 相同 ， 


E 成 权重 向 量 ， 首先 , 均匀 分 布 的 
标 空间 中 搜索 方向 是 均匀 分 布 。 其次, 每 个 
权重 向 量 和 理想 点 让 每 个 个 体 趋向 Pareto 前 沿 .在 一 定 程度 上 ， 


标 函 数值 相同 ， 则 
标 函 数值 排序 。 一 直 持续 到 种 群 中 的 所 有 个 体 都 被 
它们 的 顺序 可 以 是 
任意 的 。 对 于 这 个 经 过 排序 的 种 群 P， 如果 m <n, 解 pn 将 永远 


end for 


return F; 


择 个 体 来 引导 进化 


另外 ， 本 文 的 差分 进化 算法 需要 从 NDbest 解 集中 随机 选 
过 程 ， 在 NDbest 解 


中 的 非 支配 解 集 来 引 


SURE, XL 


Bm 
n 


时 考虑 了 NDbest 解 集中 一 些 较 好 的 解 来 增 


加 种 群 的 多 样 性 。 NDbest 解 集 将 会 用 于 下 面 的 多 策略 变异 操作 。 
2.3 多 策略 变异 操作 


差分 进化 算法 中 变异 策略 会 影响 整个 种 群 的 更 新 , 因此 选 


择 合适 的 变异 策略 
异 策略 针对 种 群 的 
不 能 推进 收敛 过 程 


besU1”， 可 能 又 因为 多 样 性 差 从 而 导致 种 群 早熟 收敛 。 


是 改善 整个 算法 性 能 的 关键 “DE/rand/1” 变 
多 样 性 有 很 好 的 效果 ， 但 是 收敛 性 不 足 可 能 
， 而 如 果 选 择 收 敛 性 强 的 策略 “DE/ rand -to- 
因此 ， 


在 选择 合适 的 变异 


策略 上 ， 尽 可 能 兼顾 其 收敛 性 和 多 样 性 ， 从 


而 有 效 的 推进 进化 


的 过 程 。 


本 文 提出 了 


种 多 策略 变异 。 在 进化 初期 ， 让 DE/rand/l 


变异 策略 发 挥 其 全 


局 搜索 的 能 力 增加 种 群 的 多 样 性 ， 随 着 进化 


不 会 被 解 ps 支配 ， 因 为 pm 中 至 少 存在 一 个 目标 ， 其 值 小 于 pr 
中 相同 目标 。 表 明 这 两 种 解 之 间 只 存在 两 种 可 能 的 关系 : Du xc 


inn 
pius 


的 多 样 性 。 高 效 非 支配 排序 的 伪 代 码 如 下 : 


高 效 非 支 配 排 序 策略 (ENS) 


配 pn: 或 者 pn 和 pn 互 不 支配 。 本 文通 过 高 效 的 非 支配 排 
过 对 种 群 进行 排序 , 选 出 前 面 N/2 个 个 体 作为 NDbest fR, BE 
保留 了 非 支 配 解 集 ， 又 保留 了 一 些 较 好 的 解 来 增加 NDbest 解 


序 通 


F-empty; /*F 为 前 沿 的 集合 */ 


也 种 群 P 按照 第 1 个 目标 函数 值 进行 升序 排序 ， 
标 函 数值 相同 ， 则 按 第 二 个 目标 函数 值 排序 ， 一 直 
所 有 个 体 都 被 排序 ; 

forall P[n] e sorted P do 


x-size(F); 
kl; 
while true do 
M F[K] 的 最 后 一 个 到 第 一 个 解 与 比较 Pln] ; 
if F[k] 中 没有 解 支配 Pln] then 
把 P[n] 添加 到 F[K] ; 
break; 


else 
kt; 
if k>x then 
xtl ; 
把 P[n] 添加 到 新 的 前 沿 ; 
break; 
end if 
end if 


end while 


如 果 两 个 解 的 血 
持续 到 种 群 


H 


IL 


的 


过 程 的 进行 ， 本 文 提出 了 一 种 新 的 变异 策略 “DE/rand-to- 


NDbesW1”， 通 过 上 一 代 高 效 非 支配 排序 选择 出 NDbest MEK 


指导 进化 过 程 进 


步 增强 收敛 性 ， 从 而 达到 了 收敛 性 和 多 样 性 


的 平衡 ， 最 后 运用 


二 项 式 交 叉 和 多 项 式 变异 。“DE/rand-to- 


NDbest/1” 变 异 策略 如 下 : 


ViG=Xet rand - Copies 


其 中 : NDbest 表示 
的 下 标 ，rand 表示 


i s if (rand (0,D < CR)orCj = j,,4). j 21,2,...D 
Jg 


j 
Xic 


=x p) t Fx, eto) (11) 
当前 种 群 高 效 非 支 配 排序 最 优 集中 随机 个 体 


0 到 1 的 随机 数 。 二 项 式 交叉 如 下 : 


(2) 


otherwise 


其 中 CR 是 交叉 概率 ， ji 73 1 到 D 中 随机 产生 的 一 个 整数 ， 


以 保证 尝试 向 量 y 


"m 
yo = 
Xi 


1- (2-2* rand)" 


: pa if (rand (0,D < p, ) 


与 当前 个 体 x 不 相同 。 多 项 式 变异 公式 如 下 : 


(13) 


otherwise 


i (14) 


otherwise 


1 
, | (2+rand)"™ —1  if(rand € 0.5) 
oj = 


其 中 p 为 变异 比率 , b; 和 a 分别 是 第 i 个 决策 变量 的 上 下 界 ， 
N 是 分 布 指数 。 多 策略 变异 的 伪 代 码 如 下 : 


多 策略 变异 操作 


if rand <1- (G / G,, ^ /*rand 是 0 到 1 随机 生成 的 随机 数 */ 


Vig = X, g t F(X, 6 — X, 5) /*DE/rand/1*/ 


else 


Vig —7XQg* rand * (Xypys — X.) ag Fx, = X.) 


/*DE/rand-to-NDbest/1*/ 


end if 


根据 式 〈12) 进行 二 项 式 交 又; 


最 后 根据 式 (130 (14) 进行 多 项 式 变异 。 
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2.4 MODE/DMSM 算法 

MODE/DMSM 算法 总 体 流 程 如 下 :a) 初 始 化 种 群 ，b) 采 用 
2.2 中 方法 对 种 群 进行 高 效 的 非 支配 排序 选 出 前 N/2 个 个 体 作 
为 NDbest R; c) 采 用 2.3 中 方法 多 策略 变异 产生 子 代 y: dy 
用 2.3 中 方法 改进 的 Tehebycheff 分 解 方法 更 新 种 群 中 的 个 体 ; 
6) 采 用 2.2 中 方法 对 种 群 进行 高 效 的 非 支配 排序 选 出 前 M2 个 
个 体 作 为 VDbes! 集 ; 了 如 果 满 足 停 止 条 件 , 算 法 停止 ,输出 PF, 


否则 返回 步骤 c) 继 续 执行 .MODE/DMSM 算法 的 伪 代 码 如 下 : 


Algorithm: pseudo code of MODE/DMSM 


随机 初始 化 种 群 P， 设 置 种 群 数量 N 和 最 大 评估 次 数 MA4X_FES, 计 算 群 


体 每 个 个 体 适应 度 值 和 每 个 目标 最 优 适应 度 值 Z = (ZiZa) > BELA 


HT 


根据 文献 [18] 可 知 ，MOEADDE 的 时 间 复 杂 度 为 O(MNT) US], 
T TEN HT 5S NX MOEADDE 的 时 间 复 杂 度 为 O(MN”) 。 
MOEADDE 在 种 群 进化 中 始终 采用 单一 的 变异 算 子 ， 而 
MODE/DMSM 在 MOEADDE 基础 上 增加 了 多 策略 变异 算 子 没 
有 额外 增加 算法 时 间 复 杂 度 ; 另外 , MODE/DMSM 运用 了 高 效 
非 支配 排序 , 根据 文献 [17] 可 知 , 高 效 非 支配 排序 所 需要 的 最 好 
情况 的 时 间 复 杂 度 是 OQMN AN) ， 最 坏 情 况 的 时 间 复 杂 度 是 
O(MN?^) U7, W MODE/DMSM 的 时 间 复 杂 度 为 O(MN?) 。 


lm 


T 


3 ”仿真 实验 及 结果 分 析 


为 了 测试 本 文 算法 MODE/DMSM 的 性 能 ， 选 取 4 种 代表 
性 多 目标 算法 进行 对 比 ， 分 别 是 MOEADI9, NSGAIIIU, 


成 N 个 权重 向 量 Am CAS... AD 。 计 算 权 重 向 量 之 间 的 欧式 距离 ， 
并 且 找 出 每 个 权重 向 量 距离 最 近 的 了 个 权重 向 量 A", AT ， 其 索引 是 
B()-0,..1,). ERRANAK 6 为 0.9， 子 代 蔡 换 父 代 的 
最 大 数量 7 为 2 ，FES=N。 
通过 2.2 中 方法 高 效 的 非 支 配 排 序 CENS) 进行 非 支配 排序 选 出 N/2 个 
个 体 作为 NDbest。 

while FES<Max_FES do 


for i= 1 to N do 

随机 产生 0 到 1 的 随机 数 rand Æ 
B(i) if (rand < ó) 

{1,...N} otherwise 


成 


t 


A- 


TR 2.3 中 方法 多 策略 变异 DE 更 新 种 群生 成 新 的 解 y。 
如 果 y 超 出 了 边界 的 范围 ， 在 边界 里 随机 生成 新 值 代替 y。 


计算 个 体 的 适应 度 值 ， 更 新 Z = (Gn. 


for7= 1 to m do 
if z; > fy) then 
z -f,0) 
end if 
end for 
更 新 种 群 的 解 ,设置 c=0， 
while —c = n, | isempty(A)) do 
随机 从 4 rU ERST k; 
if /— gyl Az) S gx, | Az) 


X,—y; 


c-ctl; 
end if 
从 4 中 清除 索引 不 


end while 
FES=FES+N,; 


end for 


通过 22 中 方法 高 效 的 非 支配 排序 CBENS) 进行 非 支 配 排序 选 出 N/2 
个 个 体 作 为 NDbest 
end while 
2.5 时 间 复 杂 度 分 析 
假设 算法 的 种 群 大 小 为 Y， 目 标 数 为 M, 


领域 大 小 为 工 。 


RVEAU9I, MOEADDEU!, NSLSP!fII MOEAIGDNSU?!, 所 有 的 
算法 都 是 在 多 目标 进化 算法 开源 平台 PlatEMOP3 上 运行 。 实 验 
仿真 的 环境 都 是 在 Intel(R) Core(TM)i7-7700HQ CPU、16G 内 
存 、2.8GHz 主 频 ，win10 64 位 操作 系统 的 计算 机 上 Matlab 
R2016a 运行 。 
3.1 测试 函数 及 评价 指标 
本 文采 用 文献 [24] 提 出 ZDT1-ZDT4、ZDT624 和 文献 [25] 
提出 DTLZI-DTLZSUSff] 10 个 测试 函数 来 测试 本 文 算法 
MODE/DMSM 的 有 效 性 。 评 价 指 标 采 用 世代 距离 
GD(Generation Distance)U6! 、 反 向 世代 距离 IGD(Inverted 
Generation Distance) 7 和 超 体积 HV(Hypervolume)P8 来 评价 解 
集 的 收敛 性 和 多 样 性 。 其 计算 公式 如 下 : 


d(v, P* 

apep, p) - 22 P) m (15) 
|P] 
.d(v, P 

pne p 2s P i i (16) 
|P'| 

HV(P)- vo (Uu (2), zx. xE O, 2] (17) 


其 中 :P 为 所 得 的 Pareto 前 沿 , P* 为 最 优 的 Pareto 前 沿 , d(v, P^) 
表示 目标 空间 上 v 与 p^ 中 最 近 解 之 间 的 欧式 距离 , Vol) 记 为 
Lebesgue 测度 。 
3.2 ”实验 参数 设置 

算法 在 每 个 测试 函数 上 独立 运行 30 次 ,设置 种 群 规 模 NP 
为 100， 最 大 函数 评估 次 数 Max_FES 为 30000， 多 项 式 变异 
7=20 、p,=1/D 。 文 本 算法 MODE/DMSM 的 缩放 因子 F 为 
0.5， 交 又 概率 CR 为 0.2; MOEADDE 的 缩放 因子 下 为 0.5， 交 
又 概率 CR 73 1; 本 文 MODE/DMSM 和 MOEADDE 中 的 邻 域 
个 体 选择 概率 5 为 0.9， 子 代 蔡 换 父 代 的 最 大 数量 为 2; 
MOEAD, NSGAIII, RVEA, NSLS 和 MOEAIGDNS 中 交叉 率 
D, 7l. 7, 230; RVEA 中 惩罚 因子 变化 率 为 a=2，, 参考 点 调 
整 频率 f, =0.1。 
3.8 ”实验 结果 分 析 

MODE/DMSM 5j MOEAD,NSGAIII,RVEA, MOEADDE.、 
NSLS、MOEAIGDNS 测试 于 ZDT 和 DTLZ 测试 函数 上 分 别 在 


GD, IGD 和 HV 标准 下 
计 结 果 如 表 1-3 所 示 ， 本 文采 


的 30 次 独立 运 
Wilcoxon 秩 和 检验 ; 
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行 的 平均 值 和 标 


NIE 25 SE 


fT 


出 算法 之 间 差 异 的 显著 拆 


E, 显著 性 水 3 


“=” 分 别 表 示 其 他 算法 明显 
MODE/DMSM， 其 


中 最 好 的 结果 用 


uj 


不 同 颜 


从 表 1 中 的 数据 在 GD 标准 


MODE/DMSM 分 别 在 2 个 测试 函 


MOEAD; 


EJ 0.05, 表 中 “+”、 
明显 劣 于 和 近似 了 
色 的 字体 表示 。 
下 可 以 看 出 , 本 文 提出 
数 (ZDT6、DTLZ5) 上 近似 于 
在 2 个 测试 函数 (ZDT3、ZDT6) 上 近似 于 NSGANI; 
在 1 个 测试 函数 (ZDT6) 上 近似 于 RVEA; 在 1 个 测试 函数 (ZDT6) 


比较 突 


éé m 
= 


H 


的 算法 


上 近似 于 MOEADDE， 在 两 个 测试 函数 (ZDT3、ZDT6) 上 近似 


于 MOEAIGDNS。 在 其 余 的 测试 函 


数 下 ，MODE/DMSM 优 于 


i 


MOEAD,NSGAIII, RVEA, MOEADDE, NSLS, MOEAIGDNS 。 


Ak E, MODE/DMSM 在 8 个 测 
DTLZ5) 上 得 到 最 好 的 结果 。 


试 函 


从 表 2 中 的 数据 在 IGD 标准 


下 可 


在 8 个 测试 函 


数 (ZDT1~4、DTLZ1~3、 


以 看 出 ，MODE/DMSM 
Ji CZDT1-4, ZDT6. DTLZ1、 DTLZ 3, DTLZ5) 


上 优 于 MOEAD; MODE/DMSM 在 5 个 测试 函数 CZDTI. ZDT 


3~4、DTLZ3~4) 上 优 于 NSGAII; 


DTLZ5 


EET 


MODE/DMSM 在 7 


个 测试 


函数 CZDTI-4. ZDT6. DTLZ3. DTLZ 5) 上 优 于 RVEA; 
MODE/DMSM 在 8 个 测试 函数 CZDT1-4. DTLZ1-4) ERF 
MOEADDE; MODE/DMSM 在 8 个 测试 函数 (ZDT1~4、ZDT6、 
DTLZ1、DTLZ3~4) 上 优 于 NSLS, MODE/DMSM 在 6 个 测试 
函数 (ZDT1~2、ZDT4、ZDT6、DTLZ3~4) 上 优 于 MOEAIGDNS。 
总 体 上 ，MODE/DMSM 在 5 个 测试 函数 (ZDT1~4、DTLZ3) 上 
得 到 最 好 的 结果 。 
AR 3 中 的 数据 在 HV 标准 下 可 以 看 出 ,MODE/DMSM 
在 8 个 测试 函数 CZDTI-4. ZDT6. DTLZI , DTLZ 3, DTLZ5) 
上 优 于 MOEAD; MODE/DMSM 在 7 个 测试 函数 CZDT1-4. 
DTLZ1、DTLZ3、DTLZ4) 上 优 于 NSGAII; MODE/DMSM 在 
8 个 测试 函数 (ZDT1~4、ZDT6、DTLZ1、DTLZ3、DTLZ5) 上 优 
于 RVEA;MODE/DMSM 在 8 个 测试 函数 (ZDT1~4、DTLZ1~4) 
EF MOEADDE; MODE/DMSM 在 9 个 测试 函数 (ZDTI1~4、 


ZDT6、DTLZ1~4) 上 优 于 NSLS 和 MOEAIGDNS。 总 体 上 ， 


MODE/DMSM 在 6 个 测试 函数 (ZDT1~4、DTLZ1、DTLZ3) 上 
得 到 最 好 的 结果 


de 1 测试 函数 在 GD 标准 下 实验 统计 结果 


MOEAD NSGAIII 


4.1652e-4 9.2005e-5 


(1.19e-4)- | (3.45e-5)- 
9.4598e-4 6.3058e-5 
(6.22e-4) - | (2.31e-5)- 
1.5232e-3 6.6748e-5 
(1.94e-3)- | (5.54e-5) = 
1.7281e-3 2.9247e-4 
(104e-3)- | (1.72e-4) - 
8.1940e-4 6.2187e-5 
(2.22e-4) = | (6.02e-5)= 
3.3229e-4 2.6684e-4 
(1.46e-4) - | (1.14e-4)- 
5.0524e-4 5.0326e-4 
(2.64e-6) - | (5.75e-6) - 
1.7300e-1 3.3643e-1 
(3.98e-1)- | (8.69e-1)- 
3.7487e-4 4.9923e-4 
(1.8le-4)= | (6.50e-5)- 
1.6000e-5 2.5470e-4 


G.16e-5)- | (1.15e-4) - 


3.2302e-3 


(8.28e-4) - 


3.3958e-3 


(1.44e-3) - 


3.3952e-3 


(1.25e-3) - 


2.1164e-3 


(8.00e-4) - 


3.8055e-3 


(1.37e-3) = 


2.8930e-4 


(1.34e-4) - 


5.1379e-4 


(2.22e-5) - 


1.8290e-1 


(1.84e-1) - 


5.1017e-4 


(6.98e-6) - 


4.2376e-2 


(1.542) - 


MOEAIGDNS 


MOEADDE 


1.0803e-3 5.1423e«0 7.3084e-5 1.2125e-5 


(6.49e-4) - (9.68e-1) - (2.47e-5) - (9.29e-6) 


9.5758e-4 6.8252e«0 4.6044e-5 7.8375e-6 


(5.12e-4)- | (137e40)- (2.01e-5) - (5.78e-6) 


2.1661e-3 6.0136e«0 5.0231e-5 4.9466e-5 


(2.21e-3)- | (126e40)- (1.33e-5) = (1.14e-5) 


2.8039e-2 3.6696e-1 3.7083e-4 4.4888e-5 


(1.71e-2) - (9.17e-2) - Q.27e-4) - (1.59e-5) 


2.4356e-3 7.5549e-2 2.0133e-4 1.6643e-2 


(1.09e-2)= | (6.84e-2) - (2.96e-4) = (2.24e-2) 


2.8448e-2 5.1296e-2 4.3197e-4 1.8284e-4 


(1.28e-1) - (3.14e-2) - (1.95e-4) - (8.07e-6) 


8.1839e-4 1.2919e-3 6.8633e-4 4.9107e-4 


(6.15e-5) - (1.92e-4) - (4.44e-5) - (1.71e-5) 


4.4169e-1 5.6068e+0 1.8665e-1 2.6823e-2 


(8.35e-1) - (1.01e+0) - (5.67e-1) - (6.84e-2) 


8.6356e-4 5.2988e-3 6.2166e-4 4.7772e-4 


(4.21e-4) - (4.00e-3) - (1.62e-4) - (1.59e-5) 


8.8713e-5 2.6775e-4 3.7241e-5 6.5693e-6 


(2.49e-5) - (5.77e-5) - (1.11e-5) - (1.58e-6) 


表 2 测试 函数 在 IGD 标准 


FH 


准 下 实验 统计 结果 


Problem 


MOEAD NSGAIII 


RVEA 


MODE/ 


MOEADDE NSLS MOEAIGDNS 


DMSM 
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1.6602e-2 5.3384e-3 3.2468e-2 1.1745e-2 7.1899e+0 6.4159e-3 3.9653e-3 
ZDTI 2 
(2.12e-2)- | (540e-3- | (2.01e-2)- (5.73e-3)- | (2.06e+0) - (8.65e-3) - (3.96e-5) 
4.8749e-2 4.9106e-3 4.4291e-2 9.6554e-3 8.5175e40 1.6885e-2 3.9602e-3 
ZDT2 2 
(715e22)- | G.15e-3)= | (5.38e2)- (3.74e-3)- | (1.79040) - (2.58e-2) - (4.39e-5) 
2.8498e-2 6.0005e-2 8.2320e-2 2.8014e-2 7.6289e+0 3.6516e-2 1.1575e-2 
ZDT3 2 
(2.75e-2) - | (7.59e-2) - | (5.14e2)- (2.04e-2) - | (1.93e40)- (3.76e-2) = (3.71e-4) 
2.1374e-2 1.0630e-2 4.6528e-2 2.1136e-1 8.3338e-1 9.0156e-3 4.1580e-3 
ZDT4 3 
(126e-2)- | (2.06e-2)- | (493e2)- (1.14e-1) - (1.97e-1) - (1.43e-2) - (1.29e-4) 
7.3477e-3 3.1696e-3 2.7225e-2 3.1064e-3 6.5237e-3 3.5963e-3 3.2963e-3 
ZDT6 2 
(L46e-3)- | (2.56e-4)+ | (7.75e-3)- | (9.11e-6)+ | (1.77e-3)- (3.88e-4) - (9.80e-5) 
2.0988e-2 2.0825e-2 2.1209e-2 4.9475e-2 2.8327e-1 2.0138e-2 2.0750e-2 
DTLZ1 3 
(4.36e-4) - | (&.44e-4)7 | Q23e3)- | (655e2)- (1.96e-1) - (7.85e-4) + (7.37e-5) 
5.4467e-2 5.4497e-2 5.4508e-2 7.5988e-2 5.6350e-2 5.1570e-2 5.6600e-2 
DTLZ2 | 3 
(1.43e-6)+ | (9.15e-5)+ | (125e-4) | (830e-4)- | (7.54e-4)= (5.61e-4) + (7.36e-4) 
1.2828e+0 | 5.6226e-1 8.7055e-1 22512e40 1.7435e«1 2.7297e-1 2.5802e- 
DTLZ3 | 3 
(2.71e+0)- | ($.11e-)- | (9.05e-1)- | (4.18e+0)- | (3.85e+0)- (4.12e-1) - (5.05e-1) 
3.8164e-1 1.6816e-1 5.4487e-2 1.5330e-1 2.9700e-1 1.6585e-1 6.0017e-2 
DTLZA | 3 
G4le-1)& | Q.10e-)- | (4.54e-5)+ | (8.37e-2)- (6.82e-2) - Q.10e-1) - (1.63e-3) 
3.3778e-2 1.3344e-2 7.9796e-2 14375e-2 5.1314e-3 4.1794e-3 22831e-2 
DTLZ5 | 3 
(692e-5)- | (1.68e-3)+ | (1.26e-2)- | (9.07e-5)+ | (2.40e-4) + (5.24e-5) + (8.29e-5) 
+/-/= 1/8/1 3/5/2 2nn 2/8/0 1/8/1 3/6/1 
JE se 尖 在 =E TK srt ills 25 2-]- zc 
d 3 测试 函数 在 HV 标准 下 实验 统计 结果 
MODE/ 
Problem | M MOEAD NSGAIII RVEA MOEADDE NSLS MOEAIGDNS 
DMSM 
8.5700e-1 8.6904e-1 8.2655e-1 8.5580e-1 0.000060 8.6834e-1 8.7143e-1 
ZDTI 2 
(1.60e-2) - | (476e-3)- | (1.59e-2) - (9.74e-3)- | (0.00e+0) - (7.21e-3) - (9.01e-5) 
4.7374e-1 5.3484e-1 4.7327e-1 5.2297e-1 0.0000e+0 5.1430e-1 5.3829e-1 
ZDT2 2 
(7.90e-2)- | (7.60e-3)- | (5.53e-2)- (825e-3)- | (0.00e+0) - (3.73e-2) - (5.78e-5) 
9.8476e-1 9.5422e-1 9.0636e-1 9.7045e-1 0.000040 9.8196e-1 1.020140 
ZDT3 2 
(3.82e-2) - | (9.67e-2)- | (5.99e-2) - (3.85e-2)- | (0.00e+0) - (4.61e-2) - (1.83e-4) 
8.4079e-1 8.6296e-1 82191e-1 5.7036e-1 4.1779e-2 8.6339e-1 8.7086e-1 
ZDT4 2 
(1.60e-2)- | (1.58e-2)- | (Q.77e2)- (1.52e-1) - (4.70e-2) - (1.07e-2) - (2.57e-4) 
4.2470e-1 4.3277e-1 3.9550e-1 4.3342e-1 4.3003e-1 4.3119e-1 4.3297e-1 
ZDT6 2 
(245e-3)- | (794e-4)- | (1.16e-2)- | Q.18e-4) | (1.69e-3)- (9.93e-4) - (4.24e-4) 
1.3947e-1 1.3965e-1 1.3953e-1 1.2506e-1 4.9144e-2 1.3676e-1 1.4000e-1 
DTLZI 3 
(4.29e-4)- | (3.23e-4)- | (8.34e-4)- (2.44e-2) - (5.26e-2) - (8.88e-4) - (1.80e-5) 
7.4465e-1 7.4463e-1 7.4446e-1 7.0005e-1 7.3781e-1 7.2068e-1 7.4425e-1 
DTLZ2 | 3 
(4.54e-5)+ | (9.42e-5) | (2.90e-4) | (1.84e-3)- (1.87e-3) - (4.67e-3) - (2.98e-4) 
32185e-1 3.6040e-1 2.3777e-1 3.3787e-1 0.0000e+0 5.1204e-1 6.1475e-1 
DTLZ3 | 3 
(3.07e-1)- | Q91e-D- | (2.55e-1)- (3.03e-1)- | (0.00e+0) - (2.44e-1) - (2.71e-1) 
5.3426e-1 6.7613e-1 7.4451e-1 6.7798e-1 5.3580e-1 6.6337e-1 7.4378e-1 
DTLZ4 | 3 
(2.3le-0D)= | (126e-1)- | (L18e-4)4 | (3.63e-2)- (8.39e-2) - (1.14e-1) - (4.70e-4) 
1.2107e-1 1.2860e-1 9.9884e-2 1.2934e-1 1.3247e-1 1.3228e-1 1.2675e-1 
DTLZ5 | 3 
(2.07e-5)- | (8.25e-4)+ | (5.94e-3)- | (5.40e-5)+ | (1.62e-4)+ (6.28e-4) + (2.11e-5) 
M 1/8/1 7/1 2/8/0 2/8/0 1/9/0 1/9/0 
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ZDT、DTLZ 系列 测试 
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， 从 图 2(a)~G) 可 以 看 


测试 函数 上 能 获得 所 得 至 


PF 
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pui 


求解 ZDT 和 DTLZ 问题 Pareto 前 沿 分 布 图 来 
直观 的 了 解 改 算法 的 性 能 。 


图 2 展示 了 MODE/DMSM 算法 在 
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[分 布 均匀 。 


上 某 次 独立 运行 的 Pareto 前 沿 分 布 


H MODE/DMSM 在 ZDT、DTLZ 系列 
的 非 支 配 解 集 能 够 较 好 的 逼近 真实 的 
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(a) ZDT1 上 的 仿真 实验 


MODE/DMSM on ZDT3 
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(c) ZDT3 上 的 仿真 实验 


MODE/DMSM on ZDT6 
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(e) ZDT6 上 的 仿真 实验 


MODE/DMSM on DTLZ2 


(gDTLZ2 上 的 仿真 实验 


算法 无 论 从 GD IGD 还 是 HV PR 
RVEA, MOEADDE, NSLS, MOEAIGDNS 的 收敛 性 和 多 样 
要 好 。 MODE/DMSM 算法 在 处 理 多 目标 问题 有 较 好 的 性 能 , 进 
一 步 体现 了 算法 的 有 效 性 。 
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终 上 所 述 ,针对 ZDT fI DTLZ 测试 函数 上 ,MODE/DMSM 
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(bZDT2 上 的 仿真 实验 
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MODE/DMSM on DTLZ4 
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G)DTLZA 上 的 仿真 实验 
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针对 多 目标 问题 , 本 文 提 出 了 一 种 基于 分 解 和 多 策略 变异 
的 多 目标 差分 进化 算法 MODE/DMSM。 本 文 运用 一 种 改进 
Tchebycheff 分 解 方法 将 一 个 MOP 分 解 为 一 组 单 目标 优化 子 问 
题 ， 通 过 标准 边界 交叉 (NBD 技术 生成 一 系列 均匀 分 布 的 权 习 


pin 


并 采用 了 高 效 非 支配 排序 有 效 保 存 了 非 支配 解 和 一 部 分 
较 优 的 解 ， 达 到 增加 收敛 性 同时 增加 其 多 样 性 ， 同 时 ， 所 提出 
多 策略 在 进化 前 期 有 效 保持 了 种 群 的 多 样 性 ， 后 期 达到 了 收敛 
性 和 多 样 性 的 平衡 。MODE/DMSM 在 多 目标 测试 函数 集 ZDT 


BÆ, 


fil DTLZ 的 10 个 测试 函数 上 进行 了 仿真 实验 ， 并 与 MOEAD、 
NSGAIII, RVEA, MOEADDE, NSLS 和 MOEAIGDNS 算法 进 
行 了 对 比 ， 结 果 表 明 所 提出 的 算法 在 求解 两 个 或 三 个 目标 的 多 
目标 问题 时 在 收敛 性 和 多 样 性 具有 一 定 的 优势 。 
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